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A SPECTRAL REFLECTANCE ESTIMATION TECHNIQUE 
USING MULTISPEGTRAL DATA FROM THE 
VIKING LANDER CAMERA 

Stephen K. Park and Friedrich O. HucL 
Langley Research Center 

SUMMARY 

A technique is formulated for constructing spectral reflectance curve estimates 
from multispectral data obtained with the Viking lander camera. The multispectral data 
are limited to six spectral channels in the wavelength range from 0.4 to 1.1 pm, and most 
of these channels exhibit appreciable out-of-band response. The output of each channel is 
expressed as a linear (integral) function of the (known) solar irradiance, atmospheric 
transmittance, and camera spectral responsivity and the (unknown) spectral reflectance. 
This produces six equations which are used to determine the coefficients in a represen- 
tation of the spectral reflectance as a linear comb. nation of known basis functions. The 
natural cubic spline basis functions generally produce the best estimates. 

Reflectance estimates are produced for a varietj' of materials that can be reasonably 
expected to occur on Mars. In each case the dominant reflectance features are accurately 
reproduced, but small period features are lost due to the undersamplmg inherent with the 
limited number of channels. The results also suggest that the estimation techhique may 
be a valuable aid in selecting the number of spectral channels and their responsivity 
shapes when designing a multispectral imaging system. 

INTRODUCTION 

The Viking lander cameras (refs. 1 and 2) will return multispectral images of the 
Martian surface with resolution higher by four orders of magnitude than has been pre- 
viously obtained. It is desired to extract spectral reflectance curves from these data. 
However, the data are limited to six spectral channels, and most of these channels exhibit 
out-of-band response. 

It is inappropriate to generate a data point for each channel by associating a reflec- 
tance value with a distinct wavelength; this is particularly true for those chaimels with 
appreciable out-of-band response. It is unlikely that data points so constructed will lie 
on the true reflectance curve and that any method of fitting a curve to these points will 
adequately approximate the true reflectance. 


Instead, the output of each channel can be expressed as a linear (integral) function 
of the unknown spectral reflectance and of the known solar irradiance, atmospheric trans- 
mittance, and catpera responsivity. This produces six equations which can be used to 
determine the coefficients in a representation of the spectral reflectance as a linear com- 
bination of known basis functions. Several sets of basis functions are discussed - poly- 
nomials and cubic splines - and their appropriateness is evaluated by constructing reflec- 
tance estimates for a variety of materials that can be reasonably expected to occur on the 
Martian surface and for the reference test charts that are located on each lander. The 
quality of these estimates is determined by comparison with the estimates that could be 
obtained by using an idealized camera response, represented by impulse functions, which 
provides an accurate reflectance sample at a distinct wavelength for each channel. 

SYMBOLS 


A 

A-1 

A 

aij 

bi 

b 

C(\) 

C"(k) 

*^i 

F(X) 

fi(^) 

hjW 


6X6 system matrix 
inverse of A 

8X8 matrix formed from A (see table n) 
elements of matrix A 
spectral sample for ith channel 

column vector with elements 0, b^, , bg, . . ., bg, 0 

cubic spline basis function 

second derivative of C(\) 

system calibration constant for .ith channel, V 

measure of error (see eq. (10c)) 

relative standard deviation of reflectance estimate 

system characteristic function 

basis function 
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h(X) 

n 


S(X) 

Ti(X) 

ti 

Vi 


X 


A 

6(X-Xi) 


eix) 


X 

Xi 

pix) 

p^(x) 

<p(x» 


column vector witli elements IiqW, Ii 7 (x) 

number of basis functions 
channel responsivity, A/W 
solar irradiance, kW/m^-^nj. 

-1 

system transfer function for ith channel, pm 
scaling constant (see eqs. (5)) 
channel voltage, V 

channel voltage corresponding to reference reflectance, V 
coefficient of basis function h^{X) 

column vector with elements sjg, x^, X 2 , . . Xg, x,^ 

spacing between adjacent spline knots (0.12 pm) 

delta (impulse) function located at Xj, 

estimate error (see eq. (8)) 

incidence angle, deg 

wavelength, pm 

location of impulse system transfer functions, pm 
location of cubic spline knots, pm 
actual spectral reflectance 

actual spectral reflectance of reference test-chart patches 
estimated spectral reflectance 


standard deviation 


ff 

(Tp standard deviation of reflectance estimate 

atmospheric transmittance 

ci 

T (a) camera optical throughput 

( ) transpose of ( ) 

CAMERA SYSTEM 
Description 

Figure 1 shows a cutaway view of the Viking lander camera. Light that passes 
from the scene through the window is reflected by a scanning mirror to an objective lens 
which images this light onto a plane containing an array of photosensors, each covered by 
a small aperture that defines the instantaneous field of view. The photosensors convert 
the light into an electrical signal which is amplified, sampled, and quantized for digital 
transmission. The camera rotates a small step in azimuth between elevation line scans 
until the entire scene of interest is scanned. General design and performance character- 
istics are listed in table I. 

The camera features an array of 12 silicon photodiodes, consisting of 4 broadband 
channels with selectable focus for high-resolution imaging, 1 broadband channel for low- 
resolution surveys, 6 narrowband channels for multispectral color and near-infrared 
imaging, and 1 narrowband red channel for scanning the Sun. Typical responsivities for 
the six narrowband channels are plotted in figure 2. The three spectral channels in the 
visual (0.4 to 0,7 pm) range were selected to produce conventional color images, and the 
three channels in the near-infrared (0.7 to 1.1 pm) range were selected to detect reflec- 
tance minima near 0.95 pm due to Fe2+ electronic transition absorption bands, 

A multispectral (color or IR) imaging mode is effected as follows : Three spectral 
channels (either blue-green-red or IR1-IR2-IR3) are alternately selected with successive 
elevation scans while the azimuth servo is inhibited. After the third elevation scan is 
completed, the azimuth servo rotates the camera a small increment and the same three 
spectral channels are again line scanned. An azimuth sweep in this manner produces a 
three-channel image of the scene. A second azimuth sweep of the scene with the other 
three channels is required to produce a six-channel multispectral image. 

The responsivity shapes of the six spectral channels depart significantly from the 
shapes that were hoped for. Photosensor-array size constraints and heat- sterilization 
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requirements dictated that interference-type filters be used. These filters have multiple- 
order transmittance peaks and are most readily made with very narrow transmittance 
bands (about 0.01 p.m for visible wavelength filters). However, the Viking camera filters 
were desired to have transmittance half -widths from 0,05 to 0.1 Mm and an average pass- 
band transmittance of over 65 percent to provide adequate sensitivity. Blocking of the 
multiple -order transmittance peaks over the silicon-photodiode responsivity range from 
0.4 to 1.1 Mm was only partially successful, and, consequently, some of the channels exhibit 
severe out-of-band leakage (fig. 2) and a generally undesirable responsivity shape. 

The undesired signal components due to out-of-band responsivity vary apnreciably 
between channels and may further increase for some of the near -infrared channels by the 
time that the landers have reached Mars. The cameras are exposed to a steady low flux 
of neutrons and gamma rays from the lander radioisotope -thermoelectric generators dur- 
ing the nearly 1-year cruise from Earth to Mars. Radiation damage by the neutrons 
degrades the photodiode performance, particularly in the infrared region from 0,8 to 
1.1 Mm where the responsivity could be reduced by an estimated 15, to 40 percent. 

Figure 3 shows a simplified circuit diagram of the photosensor array (PSA) and 
video processing electronics {VPE), Each one of the photodiodes (except the Sun diode) 
has its own preamplifier to convert the small output current into a usable signal. The out- 
puts from the preamplifiers are summed into a buffer amplifier, but only one preamplifier 
is turned on by the camera logic at any time. The gain of the buffer amplifier, or channel 
gain, is determined by the input resistances. These resistances compensate for different- 
photodiode aperture sizes (i.e., instantaneous fields of view) and filter transmittances so 
that the PSA output remains approximately constant (about 1.46 V) for all chanhels when 
the camera views the average Mars radiance. (See ref. 2.) 

Analytical Model 

The reflectance of a diffuse surface is a function of wavelength X and of the lighting 
and viewing geometry. If the blue-green-red and IR1-IR2-IR3 images are obtained at 
nearly the same Sun angle, the lighting and viewing geometry is fixed and the multispectral 
image data can be used to estimate the surface reflectance as a function of wavelength. 

That is, the dependence on lighting and viewing geometry can be ignored and the surface 
radiance N(x) can be expressed as 

N(X) = ls(X)7^(X)p(X) (1) 

where S(x) is the solar irradiance, t_(X) is the atmospheric transmittance, and p(X) 

o, 

is the surface reflectance which is, for each wavelength, the ratio of reflected light to 
incident light. 
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The camera photosensor-array output voltage due to a spatially uniform radi- 
ance input is ;i ef. 2) 

'"l = {^f y” N w W RjW dX (2) 

where D is the lens aperture diameter, /3- is the instantaneous field of view, R-. is 

^ i.1. 

the preamplifier feedback resistance, G. is the channel gain, k„. is a radiometric cal- 
ibration constant, is the camera optical throughput, R^(X) is the responsivity 

(fig. 2), and the subscript i denotes the channel -dependent parameters (i = 1, 2, . , ,, 6), 
It is convenient to combine equations (1) and (2) and express equation (2) in the simplified 
form 




1 . 2 , . 



The parameter c^ is a channel-dependent system calibration constant given by 


( 3 ) 




(i = 1, 2, . . ., 6) 


( 4 ) 


and T.(?l) is a channel -dependent system transfer function given by 


T.(\) = 


Six) tJx) r(x) R.(A) 


t, 


(i = 1, 2, , . 6) 


where tj^ is a scaling constant chosen so that 

r T,(x) dx = 1 


(i = 1, 2, . . ., 6) 


(5 a) 


(5b) 


Figure 4 presents a plot of typical values of S(X), t^(X), and t^(X). These plots 
combined with the photosensor-array responsivities plotted in figure 2 yield the six 
channel system transfer functions T^(X) plotted in figure 5. 


Calibration 


It is possible to estimate the calibration constants Cg, • • Cg by using 

equation (4); however, at least one factor — the radiometric calibration constant k - —is 
expected to change from its prelaunch calibration value by the time that the lander reaches 
the Martian surface. Thus, it is preferable to determine each constant Cj^ by using 
multispectral image data of a reference test chart (RTC). An RTC is show in figure 6; 
it provides 11 gray patches and 3 color patches (as well as 3 tribars to check the camera 
frequency response). 
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The reflectance at normal light incidence for all 14 RTC patches is shown in fig- 
ure 7. These reference reflectances (denoted by pj.{X)) have a 2-sigma accuracy of 
±8 percent. Within this accuracy the reflectance at nonnormal incidence is approximately 
Lambertian and is given by p^{x) cos l for most incidence angles (20° S t. g 60°). The 
reflectances of the gray patches relative to each other are accurate to within ±4 percent 
for any incidence angle and wavelength. 

The RTC patch reflectances p^{X) with corresponding voltages V^, Vg, . . Vg 
that fall within the selected camera dynamic range satisfy equation (3); that is, 

vf = c.(cos t) \ T.{X) pj\) dX (i = 1, 2, . . 6) (6) 

Equation (6) provides a means of estimating each calibration constant c-. That is, 
is the slope of the least-squares linear fit to the data points generated by the RTC patch 

reflectances and plotted as (cos l) \ T.(X) p^(x) dX against V[. Any significant depar- 

J 0 

ture of the data from the least-squares line gives warning that the RTC patch reflectances 
may have been changed by dust adhesion or abrasion, 

THE REFLECTANCE ESTIMATE AS A LINEAR COMBINATION 
OF BASIS FUNCTIONS 

Formulation 

Equations (3) describe the relationship between the six voltages Vgj . . Vg^ 

and the unknown reflectance p(X) of a material on the Martian surface. This section is 
concerned with formulating a technique whereby these six voltages can be used to estimate 
p(X). 

Let hj^(X), h 2 (X), . . h^(x) be a set of n known (linearly independent) basis 

functions with each defined for 0iX<«, andlet x^, Xg, , . x^^ be a set of (con- 

stant) coefficients. The function 

n ‘ ^ 

(p(x)) ^ h.(X) (7) 

j=l 

is an estimate of p(X) provided the coefficients are chosen to minimize the error 

e(x) = p(X) - <p(x)) . / (8) 
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If both sides of equation (8) are multiplied by T.(x) and integrated, ‘he following set of 
six eg- .ations (one for each channel) are obtained 

n 

= t>i - ^ ai .X. (i= 1, 2, . . 6) (9) 

j-1 

where 

V, p«s 

and 

e. = f T.(X)e(X)dX (10c) 

1 J 0 1 

The error terms provide a quantitative measure of the accuracy with which the esti- 
mate (p(X)) represents the actual spectral reflectance p(X). The coefficients 
Xj^, Xg, . . should be chosen to minimize each term je^^j. 

In general, the best estimate (p(x)> of the actual spectral reflectance /i(X) is 
obtained if e. = 0 for i = 1, 2, . . 6. Mathematically, this means that both p(X) 

and (p(X)) satisfy equations (3); and, physically, this means that two materials with 
reflectances p(X) and (p(X)> generate indistinguishable camera multispectrai 
responses. Setting e^^ = 0 for 1=1, 2,. . 6 in equations (9) yields the following set 

of six linear equations with n unknowns Xj^, Xg, . • 

n." ^ ■■ „ 

t>i=^ ai.x, (i - 1, 2, . . ., 6) (11) 

3=1 " 

If n = 6 (i.e., the number of basis functions is equal to the number of spectral 
channels) and if the resulvant 6x6 matrix A with elements a. . is nonsingular (i.e,, 
the inverse matrix A exists), then the set of equations (11) has a unique solution 
Xj^, Xg . . Xg and the (indistinguishable) estirnate |p(X)) is uniquely determined by 
equation (7). The requirement that the matrix A be nonsin^Tar is not restrictive. It 
simply means that the basis functions hj^(x), hg(x), . . ,, hg(X) and the system transfer 
functions Tj^(X), T^(\), . , Tg(x) must be chosen in such a fashion that each reflec- 
tance p(x) has one and only one ihdisttnguishable estimate (p(X)), 

8 - ■ ■■ 


Basis Functions 


The choice of basis functions is critical. Three appropriate choices are discussed 
in this section. In each case estimates (p(X)) are constructed for eight materials that 
can be reasonably expected on the Martian surface as well as for the three RTC color 
patches and one RTC gray patch (No. 5). The true (normal) reflectance for these eight 
materials is presented in figure 8, 

Polynomials. - If the basis functions are chosen to be h|(X) = 1, h 2 (x) “ X, 

. . ., hg(x) = X^, then the estimate 

6 

<P(X)> = ^ x.X^'^ 
j=l 

is a fifth degree polynomial in the variable X with coefficients xj. The corresponding 
elements a- of the square matrix A are given by 

poo . - 

a.. = T.(X) X^"^ dX 

The matrix is nonsingular. 

Figures 9 and 10 present the indistinguishable fifth degree polynomial estimates 
for the reflectance curves of the eight test materials and the four RTC patches, respec- 
tively. In each case the estimate is plotted for the effective range of the six channels, 
namely, 0.4 ^ X ^1,1 pm. The estimate is show as a continuous curve and the actual 
reflectance as a sequence of smaU circles. 

For computational purposes the 6 system tranofer functions and 12 test reflectances 
were measured at 71 equally spaced wavelengths X = 0.4, 0,41, 0.42, . . ., 1.09, 1.1. 

Each system transfer function Ws assumed to be zero Outside the interval 0.4 i X s 1.1, 
and, consequently, all integrals were limited to this range, Simpson's rule was used 
throughout to evaluate the coefficients a^^^ given by equation (lOa) and the spectral sam- 
ples b-^ given by equation (10b), 

The polynomial estimates generally agrevs well with the actual reflectances in a 
reduced domain, say 0,50 ^ X fi 1.00. Least agreement occurs when the reflectance is 
most complex, that is, for the three RTC color patches. There is a pronounced tendency 
for the estimation error je(X) | (eq. (8)) to become large as the endpoints (X = 0.4 
and X = 1.1) are approached. 

Cubic splines .- A more modern and usually better method for fitting curves to data 
involves the use of a cubic spline. A cubic spline is a curve defined on a sequence of sub - 
intervals; the endpoints of these subintervals are known as knots^ On each subinteryal 
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the curve is a cubic polynomial, and at the knots the adjacent polynomials are joined in 
such a fashion that the resulting curve (the cubic spline) is continuous with continuous 
first and second derivatives. 

The knots used for all the calculations in this paper are located at the wavelengths 
\j = 0.45 + (j - 1)A (j = 1, 2, . . 6) 

where t?:ie spacing is A = 0.12 pm, The basis functions are given by ^ref, 3) 
hj(X) = (j = 1, 2, . . ., 6) 

where C(A.) is the function defined by (see fig. 11(a)) 


6A^ L 


A^ + 3A V - IM) + 3A(A - Ul)^ - 3(A - 1^1)^ (|X| 5 A) 


C(^) =< 


:(2A .. !X!)^ 


6A" 


(a < Ixl < 2A) 


(IXI S 2A) 


On each subinterval (of width A), CfX-X.) is a cubic polynomial in X; and at 
\ V J/ 

X-X-) is continuous with continuous first and second derivatives. Con- 
jy 

sequently, the estimate ^ Xj C^X~Xj^ also has these properties and so it is a cubic 

j=l 

spline. The six basis functions are shown in figure 11(b), The elements a- of the 
square matrix A are given by 




( 12 ) 


The matrix is nonsingular. 

Figures 12 and 13 present the indistinguishable cubic spline estimate for each of the 
eight materials and four RTG patches, respectively. The estimates agree well with the 
actual reflectances on a reduced domain, say 0.45 ^X 5 1.05. As with the polynomial 
estimates, the least agreement occurs for those reflectances which have the most complex 
shape. There is still a tendency (but less pronounced) for | e(?i.) | to become large as the 
endpoints are approached. 
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Natural cubic splines. ~ The cubic spline reflectance estimates can be improved by 
adding basis functions. Specifically, figure 11(b) reveals that, at any point in the interv,tU" 
0.45 g A. g 0.57 and 0.93 g X g 1.05, only three basis functions contribute to the estimate; 
whereas four basis functions contribute at any point in the interval 0,57 g A g 0.93. This 
situation can be corrected by adding two basis functions hQ(A) and h,^(A) centered at 
Aq = 0.33 and A,^ - 1.17, respectively. The resulting estimate becomes 

7 

<p(A)> = ^ X. c(a-A.) (13) 

3=0 

The addition of the coefficients Xq and x,^. modifies equation (11) to be 
7 

bf = ^ ^ij^j (i = 1, 2, . . ,, 6) (14) 

3=0 

where the coef/^-vii'nls a^^ are given by fjuation (12). To determine the eigM coefficients 
Xq, Xj, . = .. it is necessary to add two equations to the six given by equations (14). 

A natural way to do this is to require that the cubic spline estimate given by equation (13) 
have a zero second derivative at the knots A^ = 0.45 and Aq - 1,05. The estimate then 
becomes what is known as a natural cubic spline. 

It follows from figure 11 that the two additional equations are 
Xq C”(A) + x^ C"(0) + Xg C'(-A) = 0 
and 

Xg C"(A) + Xg C"(0) + x,^ C"(-A) = 0 
which can be written, respectively, as 

Xq - 2xj^ + Xg = 0 (15a) 

and 

Xg - 2Xg + x,y = 0 (15b) 

The 8x8 matrix A formed from A by adding rows corf esponding to equations (15) is 
compiled in table H. The matrix is nonsingular. 
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If X and b denote the column vectors with elements Xq, x^, . . and 

0, bg, 0, respectively, equations (14) and (15) can be written as 

Ax = b (16) 

The indistinguisiable natural cubic spline estimate corresponding tc the samples 

b^, bg, . . bg is given by equation (13) where Xq, x^, , . are determined by 

solving equation (16). 

Figures 14 and 15 present the indistinguishable natural cubic spline estimate for 
each of the eight materials and four RTC patches, respectively. The estimates agree 
quite well with the actual reflectances on the entire domain 0.40 § \ « i. lo. The agree- 
ment is generally better than that obtained with the fifth degree polynomial and cubic 
spline estimates, particularly at the endpoints. Thus, the results indicate a preference 
for the use of natural cubic spline reflectance estimates. 

Notice specifically in figures 14 and 15 that, for the six relatively simple reflec- 
tances (i,e., pinacates 5 and 28 A, Syrtis Major, augite, average Mars, and the RTC gray 
patch), the natural cubic spline estimates are excellent. Even for the three more complex 
reflectances (i.e., limonite, hypersthene, and olivine), the estimates are very good. The 
dominant features are reproduced; however, small period features are lost. Note par- 
ticularly in figure 14 that the dominant absorption band (at X ~ 0.93 |um) for hypersthene 
is quite accurately estimated. For those three materials in figure 15 with a very complex 
reflectance (i.e., the RTC blue, green, and red patches), the shape and location of the most 
dominant features are clearly identifiable. Thus, the RTC can serve not only as an aid for 
calibrating the cameras on Mars but also as a check on the validity of the estimation pro- 
cess for multispectral data received from Mars^ 

THE REFLECTANCE ESTIMATE AS A LINEAR COMBINATION 
OF SYSTEM CHARACTERISTIC FUNCTIONS 

It is possible to derive from equation (13) an alternate, mathematically equivalent, 
but generally preferable representation for the natural cubic spline estimate. This repre- 
sentation is 

6 .. 

(p(\)) = ^ b. f.(\) (17) 

. - i= 1 

where b^ = Vj^^c^ and f j_(Ji), fgt^)? • • •} f-Q{X) are system characteristic functions. 
That is, fj^(X) is the (unique) indistinguishable natural cubic spline reflectance estimate 
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corresponding t6 the (idealized) multispectral sample with = 1 and = 0 for 
j i. These six functions are presented in figure 16. 

Images often cover spatially nonhomogeneous scenes. It is then desirable to extract 
reflectance estimates for various regions of the scene with a single set of system transfer 
functions and calibration constants Cj. The system characteristic functions are, 

therefore, fixed, and all that changes from one region of the scene to the next is the multi- 
spectral sample ^b^, bg, . . bgj. In this situation, equation (17) is the most efficient 
method to calculate reflectance estimates for these various regions. 

Although each system characteristic function is determined by all six transfer func- 
tions, each f.(x) is dependent primarily on the corresponding system transfer fuiic- 
tion T^(X). It is possible, with equation (17) and figure 16, to observe at each wavelength 
which channel (or channels) provides the primary contribution to the reflectance estimate. 

Equation (17) also facilitates the assignment of a statistical confidence interval to 
the reflectance estimate (p(\)) . Specifically, suppose the transfer functions are known 
with accuracy but the true values of bj^, bg, . . ., bg are uncertain. (Each b^ must 
be estimated from digital data that are contaminated by electronic and quantization noise*) 
Then, the system characteristic functions are known with accuracy but each is a ran- 
dom variable with mean b^ and variance By assuming also that b^ bg, . . ., bg 

are independent, the reflectance estimate is 

6 

(p(X)) = ^ b. f. (\) 
i= 1 

and the variance of the estimate % (^) is 
i=l 

2 2 2 2 

If all six variances are equal to a common value cr , then o'p W = cr F (x) where 
the relative variance of the estimate is 

6 

F^(x) = ^ f.2(x) (18) 

i=l 

At each wavelength, the statistical confidence interval centered about (p(X)> has a width 
proportional to the relative standard deviation of the estimate F(X). Figure ■.17 is a plot 
of F(X) against X; the most reliable reflectance estimates occur at those wavelengths 
where F(x) is smallest. 
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To derive equation (17), recognize that equation (IS) can be written as 
<p(X)> = x'*' h(X) 

where ( )^ denotes transpose and h(X) is the column vector with components c(x-Xq), 
c(x-Xi), . . c(x-X,y). Since A is nonsingular, the inverse A”^ exists and x can 

be written as 

X = A"^b 

Combining these last two equations yields 
<P(X)> = b’^(A'^)’^h(X) 

The 8x8 matrix (a*^) transforms the natural cubic-spline basis set into a new 
set of basis functions fQ(X), fj(X), . . ., f,^(X) where f(X) is a column vector with the 
new basis functions as its components and 

f (X) = (aT h(x) 

Since the first and last components of b are zero, fQ(X) and fr^(X) contribute nothing 
to the estimate and 

6 

(p(X)> = b*^ f(X) = ^ b^ f.(X) 
i=l 

which is equation (17). 

THE REFLECTANCE ESTIMATE FOR A SYSTEM 
WITH IMPULSE TRANSFER FUNCTIONS 

The technique discussed herein produces reflectance estimates when the system 
transfer functions have an irregular shape with appreciable out-of-band transmittance. 

It is informative to observe what this technique reduces to in the idealized situation where 
each system transfer function can be represented as an impulse function, thereby producing 
a data point for each channel by associating a reflectance value with a distinct wavelength. 

Suppose that 

Tj(X) = o(x-xJ (i = 1, 2, . . ., 6) 
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where the impulse system transfer functions occur at the discrete wavelengths 
X21 . • •) Xg. Equations (10a), (10b), and (10c) become, respectively, 

'>i = lo” pW ^ = p(h) 

and 

The elements a^^ and samples b^ are then the actual values of the basis fimctions h^(X) 
and reflectance p(X), respectively, at the wavelengths Xj. For an indistinguishable 
-reflectance estimate, e- = ^(^i) = 0 so that the reflectance estimate interpolates the 
data points p(X^), p^X2)> ■ . p(^0)* That is, if, for example, the natural cubic 

spline basis functions were used, the reflectance estimate would be a natural cubic spline 
which passes through the points p(Xj^, p(X2^, . . p(Xg). 

Figures 18 and 19 present the reflectance estimates and figures 20 and 21 present 
the system characteristic fimctions and relative standard deviation of the estimate, 
respectively, for an idealized system with impulse system transfer fimctions located at 
the six interior knots of the natural cubic spline, that is, 0.45, 0.57, 0.69, 0.81, 0.93, and 
1.05. If a six-channel system with impulse system transfer functions were designed to 
reproduce estimates best for a family of reflectances at all v^velengths 0.4 i X i 1.1 pm, 
the impulses should probably be evenly spaced and located at or near these six wavelengths. 

Comparison of figures 14 and IS with 18 and 19 reveals that both sets of estimates 
are very similar and that neither is consistently superior. This observation is supported 
by the data in table m which presents the root-mean-square (rms) error between the actual 
reflectance p(x) and the estimated reflectance (p(x)) corresponding to figures 14, 15, 
18, and 19. Three important conclusions can be made: 

(1) The proposed technique for constructing reflectance estimates completely com- 
pensates for system transfer functions with an irregular shape and appreciable out-of-band 
transmittance. 

(2) The loss of small period features in the reflectance estimate is due to the under- 
sampling Inherent in the limited number of spectral channels and is not due to the irregu- 
lar shape and appreciable out-of-band transmittance of the system transfer functions. 
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(3) System transfer functions with an irregular shape and appreciable out-of-band 
transmittance are not necessarily undesirable provided that they are properly accounted 
for by the reflectance estimation process. 

CONCLUDING REMARKS 

A teclinique was formulated for constructing spectral reflectance curve estimates 
from multispectral data obtained with the Viking lander camera. The multispectral data 
are limited to six spectral channels and most of these channels exhibit appreciable out-of- 
band response. The output of each channel was expressed as an integral function of the 
spectral reflectance producing six equations which were used to determine the coefficients 
in a, representation of the spectral reflectance as a linear combination of known basis 
functions. Several sets of basis functions were investigated. Of these the natural cubic 
spline basis functions generally produced the best reflectance estimates for a variety of 
materials that can be reasonably expected to occur on the Martian surface. 

It was found that the dominant reflectance features of these materials can be faith- 
fully reproduced; only small period features are lost due to the undersampling inherent in 
the limited number of spectral channels. The technique completely compensates for sys- 
tem transfer functions with irregular shapes and appreciable out-of-band transmittance. 
Moreover, the technique may be a valuable aid in selecting the number of spectral ch?<nnels 
and their responsivity shapes when designing a multispectral imaging system. This design 
approach should prove of value especially if spectral reflectance properties of interest are 
known a priori and if the transfer function shapes are desired to be broad to obtain good 
signal-to-noise ratios. 

The reference test charts located on the lander consist of simple gray and complex 
color spectral reflectances. It was shoylm that the shape and location of the dominant fea- 
tures for ti e complex color reflectances are closely identifiable, Thus, the reference test 
chart can serve not oiily as a means for calibrating the cameras on Mars but also as a 
means for checking the validity of the estimation process for multispectral data received 
from Mars. 

Langley Research Center 

National Aeronautics and Space Administration 

Hampton, VA 23665 . = 

September 13, 1976 
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TABLE I.- GENERAL DESIGN AND PERFORMANCE CHARACTERISTICS 
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TABLE HL- THE ROOT- MEAN- SQUARE DEVIATION BETWEEN 
THE ACTUAL REFLECTANCE AND ITS ESTIMATE 


Reflectance of - 


Average Mars . . . . , 
Syrtis Major ...... 

Pinacates 5 

Pinacates 28A 

Augite . . . . , . . . . 

Limonite 

Olivine . . . . . , . . 

Hypersthene 

RTC gray 

RTC blue . . . . , . . 
RTC green . , , . , . 
RTC red 

^See figures 14 and 15. 
^See figures 18 and l9. 


System transfer functions 
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Figure 3.- Simplified circuit diagram of photosensor array (pSA) 
and video processing electronics (VPE). 
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Figure 4.- Plots of S(A), t^(x), and t^(\) against X. 
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Figure 5.- Typical system transfer functions. 
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Figure 12.- Cubic spline reflectance estimates of the eight test materials. 
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Figure 13.- Cubic spline reflectance estimates of the four RTC patches. 
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Figure 14.- Natural cubic spline reflectance estimates of the eight test materials. 
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Figure 16.- System characteristic functions f^(X). 
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Figure 18.- Natural cubic spline reflectance estimates of the eight test materials for 
an idealized system with impulse system transfer functions. 
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Figure 19.- Natural cubic spline reflectance estimates of the four RTC patches for an 
idealized system with impulse system transfer functions. 
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Figure 21.- Relative standard deviation of the reflectance estimates for an 
idealized system with impulse system transfer functions. 
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